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POEdciom vorticity) OL C-VeCEor. 

Baogure 9. Distribution of CTD stations. 

Figure 10. A 3-D view of the stations over the shallow and 
deep waters. 

Figure 11. Density differences between pairs of stations 


of Section Db. 
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Figure 12. Density differences between pairs of deep 
stations. 

Figure 13. Linear reduction of wind stress with respect to 
depth. 

Figure 14. C,/f? values at columns 1 and 2. 

Figure 15. (C,/f* values at columns 3 and 4. 

Pigune = 16- C,/f? values at columns 5 and 6. 

Figure 17. C,/f?* values at columns 7 and 8. 

Figure 18. C,/f? values at columns 9. 

Paouises oe C,/f? values at columns (from left to tigi 
9,8 and 7. 

Figure 20. C,/f? values at sections A and B. 

Figure 21. C,/f? values at sections C and D. 

Figure 22. oy fia values at section E. 

Bcguies23, C,/f? values at sections (from left to right) 
Bh, © and A 

Figure 24. C,/f? values at surface without wind (upper 
plot) and with wind (lower plot). 

Figure 25. C,/f? values at 50 meters (upper plot) and 200 
meters (lower plot). 

Figure 26. C3/f? values at three layers (surface, 50 
meters and 200 meters). 

Figure 27. C-Vectors (for ageostrophic circulation) eae 
eolummsa9 5S cade 7 x 


Figure 28. C,/f values at Ccolummgei@andea. 
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Figure 33. C,/f values at columns (from left to right) 
Sf and 5. 


Figure 34. C,/f? values at sections A and B. 

mrogure 35. ef? values at sections C and D. 

mrgure 36. C,/f? values at section E. 

Figure 37. a values ae Sections (from left to right} 
E,C and A. 

Figure 38. W at surface. 

Figure 39. W at 50 meters (upper plot) and 200 meters 
(lower plot). 

Figure 40. Wat three layers (surface, 50 meters and 200 
meters). 

Figure 41. A satellite image of the area (15 May, 
HOS] ) . 

Figure 42. Vertical ADCP velocities at 20 meters (upper 
plot) and W function at 20 meters (lower plot). 

Figure 43. Vertical ADCP velocities at 50 meters (upper 


plot) and YW function at 50 meters (lower plot). 
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A. GULF OF THE FARALLONES NATIONAL MARINE SANCTUARY 
1. Geophysical Characteristics 
The boundaries of the Farallones shelf and slope 
region extend from Point Ano Nuevo on the south to Point Reyes 
On the north and westward from 10 km to 90 km off shore. 
This continental margin, which includes the coastal 
zone off San Francisco and is characterized by a shelf witha 


maximum width of 45 km, is the widest shelf along the entire 





Figure 1. Continental margin profile in the vicinity of 
Santa Cruz mountains (Shepard, 1941). 


west coast of the United States (Figure 1). This is, however, 
narrower than most of the shelves off the eastern and southern 
coasts of the United States. The shelf width reduces 
eeagually to the north. Another characteristic that 


distinguishes the region is that the submarine canyons are 


almost completely limited to the continental slope, beginning 
well out from the coast. 

West of Point Ano Nuevo, the direction of the 
continental slope veers northwest. The steepest part of the 
Slope branches away from the shelf border north of latitude 
37° and a relatively smooth gradient extendsmae depths of 
1,000 to or lac 0em West of the Farallones Islands, the 
steepest zone 1s close to the shelf margin, but farther north 
and beyond Cordell Bank, there is a circular shape to the 
upper limit of the continental slope. 

The fact that most canyons are restricted to the 
continental slope in this region may have some bearing on the 
relatively smooth character of the upper parts of this slope 
(Figure 2). However, there are many canyons on the outer 
slope. One of these canyons is the Pioneer Canyon between 
stations 8, 9 and 10 of section A and stations 11, 12 and 13 
of Section B. The maximum depth of this canyon is about 3,400 
m between stations 10 and 11 (Figure 3). 

According to earlier studies (Stetson, 1936 and 
Shepard et al. 1941), currents in canyons might alternate in 
direction, but more often the strongest currents are down the 
canyon. It was thought in these early studies that the period 
of alternation was not tidal. 

Shepard et al., (1979) believed that in areas with 


small tidal ranges, where short period alternations occur at 


considerable depth, the tidal current oscillations are 
PiscoGdeniatem ro us treqgional’ Current oScillations at other 
frequencies. 
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Figure 2.Profiles down the main submarine canyons off San 
Francisco (Shepard, 1941). 


An oval seamount, the Pioneer Seamount lies northwest 
of Pioneer Canyon and southwest of Farallones Islands. South 
of Pioneer Canyon is another oval seamount, Guide Seamount. 
Both rise about 1,000 m above the sea floor at the base of the 
slope, and each has more than one summit. These features may 
also play an important role in the water circulation and the 
currents of the region. 

2. Water Masses 

To study the water masses of the Farallones National 
Marine Sanctuary, the characteristic water masses in the North 
Pacific which may influence the currents west of San Francisco 


must be analyzed. 
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Longilude 
Figure 3 C.T.D. stations during the Farallones shelf and 


Slope study. 

According to Sverdrup et al. (1964) the water mass 
which dominates the upper waters of the North Pacific is the 
Subauer Te water mass. The Subarctic water mass is 
characterized as the water mass northward of latitudeusaas 
North. For the sake of convenience, this name is used for the 
entire water mass above the latitude of 239 North. This wees 


mass 1S advected eastward and then southward. 


Between 22° 2 ais Wipe iri erel depths below 300 Mm, 
temperature and salinity (T-S) curves show waters intermediate 
in character between subarctic North Pacific and the 
Beatorial water. The water at the region of interest is 
formed by the mixing of these two large and well defined water 
masses. It was determined by Collins (1992) that a poleward 
penetration of Equatorial water (associated with the under 
currents) occurred along shore over the continental slope at 
the Gulf of Farallones Islands between 16 and 21 May 1991. 
Also, the waters were found to be warm and fresh showing 
Characteristics of subarctic waters (at station 48). 

Sverdrup and Fleming (1941) and Tibby (1941) 
discovered that near the coast (30 km), the water mass had 
characteristics closer to the equatorial water mass than to 
Subarctic water mass. It is obvious that close to the coast, 
the subarctic North Pacific water mass has less effect on than 
Equatorial Water mass characteristics. 

According to Hickey (1979) the California Undercurrent 
1s a northward subsurface flow off the west coast of North 
America over the continental slope. The main flow, which 
Carries equatorial type water, was found at 200 m depth. 
During the Farallones shelf and slope study, the undercurrent 


had maximum observed velocities at a depth of 100 to 200 m. 


It has a jet-like structure, both vertically and horizontams 
and extends to the bottom over the slope. 

The most permanent current of the shelf circulation 
during upwelling is the coastal jet which is a southward flow 
with a maximum speed located between 5 km and 25 km from the 
coast. Huyer et ail. (1978) and Huyer et al. (1975) foundetiae 
the seasonally averaged jet is stronger during the spring than 
in the summer. 

According to Yoshida (1967), "processes near coastal 
boundaries respond greater and faster to time variable winds 
than those in the central oceans." During the first two weeks 
of September 1979, there were weak westerly winds which were 
not favorable for the maintenance of strong upwelling. At that 
period Breaker and Gillilard (1981) used IR images from the 
Advanced Very High Resolution Radiometer (AVHRR) to show the 
upwelled waters, near San Francisco, enter and leave San 
Francisco Bay. From the observed time sequences the upwelled 
waters appear to reflect the strong tidal influence at this 
location (Conomos, 1979). 

The Davidson current which is a northward near-shore 
flow is commonly observed north of Point Conception during the 
winter months (Reid et al., 1958), but can occur in other 
seasons under favorable wind conditions (Wyatt et al., 1972). 
The fact that drift bottles released north of Point Conception 


have been recovered as far north as Vancouver Island 


(Schwartzlose, 1963; Crowe and Schwartzlose, 1972) is strong 
evidence that the flow is usually continuous alongshore. 

The California Current (CC) which is more fully 
developed during spring and early summer has its core 
approximately 100 - 200 km offshore (Chelton, 1984) and for 
this reason, does not greatly affect the area of study which 
extends only 90 km off shore. 

3. Weather Conditions 

Weather conditions during the cruise are shown in 
Figure 4. The hourly averaged wind vectors are from west- 
northwest as usually happens during this season in this area 
(Huyer, 1983). The wind speed ranges from 3.8 to 17.2 m/sec 
and contributes to the upwelling process in the region. Air 


temperature is relatively low, between 9.89 c and 11.69 c. 


B. GENERAL CLIMATOLOGICAL PICTURE 

The data used in this thesis was gathered as part of an 
Environmental Protection Agency (EPA) funded plan to study 
oceanographic currents and hydrographic conditions in the area 
of the Gulf of the Farallones National Marine Sanctuary 
Meagure 3). 

According to the schedule, five surveys were conducted. 
This survey was carried out between May 16 to 21, 1991 aboard 
the research vessel Point Sur and was the second survey of the 
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Figure 4. Hourly averaged winds measured at 10m height on 
board R/V Point Sur during the cruise of May 16 - 21) ies 
(Jessen et al., 1992). 


During this cruise, a three-dimensional map of 
hydrographic data was obtained in order to better understand 
the currents of the area. 

Additionally, data was acquired from Advanced Very High 


Resolution Radiometer (AVHRR or AVHRR/2) instruments on board 


TIROS-N/NOAA polar orbiting satellite and interpreted as Sea 
Surface Temperature (SST) images. 
Table I. ISLANDS, ROCKS AND SHOALS IN THE GULF OF THE 


FARALLONES NATIONAL MARINE SANCTUARY. 


= 


Islands , Rocks Latitude Longitude 


and Shoals Noreen West 


1.S.E.Farallon 


3.North Farallon Is. Las 084 


5.Noonday Rock 


16.Hurst shoal 


7.Maintop Is. 





The locations of the Farallones Islands, rocks and shoals 
Mear Stations 22, 23 of Section C and 37, 38 of Section D, are 


listed in Table I. 


C. PURPOSE OF THE STUDY - CONTENTS 

This research investigates the three-dimensional 
ageostrophic pseudo-vorticity fields and the pseudo-vorticity 
fields of the total circulation near continental shelves and 
slopes extracting the necessary information from. the 
hydrographic and wind data sets (Figure 5). 

The California Current (CC), which is the eastern part of 
the subtropical North Pacific gyre and extends from Washington 
State to Baja California, is the main current at the surface. 
The California Undercurrent or Inshore Current (IC) with 
poleward flow is the main current below the surface in this 
region. The geostrophic part of these currents can be 
directly obtained from the hydrographic data sets (CTD data 
sets). 

The ageostrophic portion of these currents can be 
indirectly obtained by the C-Vector which can be calculated 
from hydrographic data and surface wind stress data. 

In this thesis, Chapter II covers the C-Vector concept and 
the techniques for the C-Vector calculations. Chapter III 
describes instruments, methods and data processing techniques. 
Chapter IV provides detailed analysis of the plots of the C- 
Vector for the ageostrophic circulation and the C-Vector for 


the total circulation (geostrophic and ageostrophic). 


ao 


CTD DATA WIND DATS 


AGEOSTROPHIC PSEUDO-VORTICLILY 
FIELDS 


PSEUDO-VORTICIFY FIELDS FOR THE 
TOVAL CIRCULATION 





Figure 5. Schematic diagram for the calculation of 
ageostrophic pseudo-vorticity fields and then the pseudo- 
femetelty fields for the total circulation. 


Finally, Chapter V summarizes the results and offers 


recommendations for future work. 
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II. C-VECTOR CONCEPT AND TECHNIQUES 


A. OBJECTIVE 

The main objective of the C-Vector concept is to Compu 
a three-dimensional pseudo-vorticity fields from the densama, 
fields and its derivatives, and then to obtain a three- 


dimensional circulations 


B. THEORY 
1. Meteorology 

The general idea was first proposed for atmospheric 
circulation with the Q-vector concept of Hoskins et al. 
(1978). Hoskins et al. found that the Q-Vector was obtained 
from the w-equation and can be explained on the synoptician’s 
chart more precisely than in earlier studies (Sutcliffe, 
1947). This theory is very helpful in diagnosing the synoptic 
and frontal vertical Gircularvon- A disadvantage of this 
concept is that it ignores the barotropic part of the 
rotational ageostrophic flow. According to Xu (1992) this 
problem can be overcome by the C-Vector concept which resolves 
the ageostrophic circulation more conveniently and allows a 
new way of solving the C-Vector equation. The C-Vector 
equation is a three-dimensional extension of the Q-Vector 


equation (QG diagnostic equation). 
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2. Physical Oceanography 

The C-Vector concept can also be applied to oceanic 
circulation (Chu, 1992). With the help of this concepewameaes 
set of hydrographic and surface wind data, the three- 
dimensional pseudo-vorticity fields ou ageostrophic 
circulation in the oceans are diagnosed. Also revealed are the 
three dimensional pseudo-vorticity fields of the total 
circulation which includes both geostrophic and ageostrophic 
flow. 

Only CTD data are needed for this calculation of the 
currents, obviating a need for more expensive and difficult 
direct measurements of ocean currents in sufficient quantity. 

There is never a perfect geostrophic balance. A 
geostrophic balance is traditionally assumed for calculation 
of the currents because the ageostrophic part is considered to 
be very small. However this assumption is of questionable 
validity especially for coastal waters of the Farallones 
National Marine Sanctuary where there are significant 
temperature and salinity gradients like the one at Section & 
and: bo (Fragure —3) = 


The thermal wind equations, 


a 9 2 -g OP , 
5 (PLY) SE Se (pfug) <9 SE (1) 
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(where p is the sea water density, f is the Coriolis parameter 
f=2Qsing) are used for the calculation of the geostrophic 
currents (Ug. Vg) . The geostrophic currents are not, however, 
the total currents. In order to have more realistic flow 
fields the ageostrophic circulation must be included. 
Following Chu (1992), the basic equations for coastal waters 
including ageostrophic circulation, wind stress and excluding 


Buoyancy are: 





a oy* 
: = = = 2 
Wren ee eV) oe (2) 
0. oy” 
ee Vee Uae) ees (3) 
(2 41-V) b+N2w=0 (4) 
Ce 
Ou, OV. OW_ 
Ee oe ee (5) 


where N is the Brunt-Vaisala frequency (assuming N 5s 
constant), f is the Coriolis parameter (on an f - plane), 
Pete v,w), and (Y’,Y") is the turbulent momentum flux due to 


), 


surface wind stress (T,,Ty 


cl 


x 
Tl pe07 ‘ leo (6) 
2] 0 
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The decomposition of the total current velocity into 


geostrophic velocity and ageostrophic circulation results in: 


V=V,+V (7) 


Using the above equation (7) and after the 
decomposition of the basic equations (2) to (4), the system of 


equations will have the form, 


H(V,) =R(V_,p, ¥", *) (8) 


where H, SR are specific differential operators. From 
functions H and ®, the three-dimensional ageostrophic 
circulation can be calculated using equation (8) together with 
CTD and surface wind data. 
3. Development of C-Vector equations (a Diagnostic 
System) 
Equations (2) to (4), for coastal waters on an £ - 
plane, can be decomposed using equation (7). The result will 


be the following equations: 
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N@w_=- (2 +V_¥) b (11) 


Cea 
Cilly Oa eel, I (12) 
Cau ao Z 


faciiGg the derivative or equation (9) with respect to 
z, and f times the derivative of equation (11) with respect to 


y, and adding them together (similarly working for equations 


(10) and (11)) results in the ageostrophic pseudo-vorticity 
equations: 
6) 0 
eae) Gi ea) =e eo) 
6) 6) 
oz (Feu) ~ a (Nfw,) =2C, (14) 
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eB ay EEE 
oa £( 2% dV, AV, dU, ee (18) 


“Ox Oy ox dy 


7, 


— 


C, = (Gaaues (19) 


The above equations yield the three components of the 
C, - Vector. Having the C, - Vector, the three dimensmaaam 
ageostrophic circulation can be founde 

Another important vector is the C-Vector (Cyr Cy, Cz) 
for the total flow, geostrophic and ageostrophic. The three 


components of this vector can be found with the following 


process. 
Using equation (13) for the total flow (VeVtVa and 
W=WgtWa but Wg=0 => wew,): 
e, re, 0 6, e, 
Se Ge — (N*w) =-—— (f° ~—— (F* (he = (20 
3, | Oe W) =~ (E*v,) 35 | V5) Sb w,) ) 
3 (21) 
tae (£°V) +2 6)=2 Cea 
ee (22) 
f-) 2 dz 
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4. Physical Significance of C-Vector 
There are two important physical meanings of C-Vector. 
a. Direction of C-Vector 
There 1S a non-dimensional formula which can 
appropriately show the relation between C-Vector and 


ageostrophic pseudo-vorticity 


VxV,=2R,C (25) 


V-V,=0 (26) 


where Ry is the Rossby number (R)=U/(fL)). 

Equation (25) shows that C-Vector is proportional 
to the ageostrophic vorticity. A C-Vector stream line can be 
considered as an ageostrophic vortex line (Figure 7). 

b. Relation between w, and vertical vorticity of C- 

Vector. 


Riis wcelatrent can be found from: 


Ws, 





Figure 7. Ageostrophic circulation induced by C-Vector in 
three-dimensional space. 


-k:Vx (equation 26) = 


Vw, =-2R,K'VxC (28) 


Equation (28) demonstrates that the vertical 


velocity (w,) 1s caused by the vertical vorticity of C-Vector 
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Figure 8. Venetealeeveloclty ow. sinduced by horizontal 
—maOicataon VOretelty) Of C-Vector. 
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Combining equation (29) with equation (28) results 


mache following: 


Vew,=-RoP (30) 

Equation (30) is the Poisson equation. If normal 
Mea-ssOlutions are considered, the vertical velocity w, has 
the same sign as WV, therefore VW will be very useful because 


its positive (negative) values show positive (negative) 


ab 


vertical ageostrophic velocities. The solutions provide 
information regarding upward (upwelling) and downward 


(downwelling) motion. 
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III. DATA COLLECTION AND PROCESSING 


The R/V Point Sur departed from the Moss Landing Harbor on 


May 16, 1991 with a combined scientific crew from Naval 


Postgraduate School (NPS), Moss Landing Marine Laboratory 
(MLML), Scientific Applications International Corporation 
(SAIC), California Academy of Science (CAS) and_e the 
Environmental Protection Agency (EPA). The purpose of the 


cruise was to study oceanographic currents and hydrographic 
conditions in the area of the Gulf of the Farallones National 
Marine Sanctuary. The ship proceeded as planned along five 
across-shore transactions 20 km apart beginning with CTD 
station 1 (Figure 3). A total of 50 CTD casts were conducted 
to within approximately 25 m of the bottom. A brief ADCP 
Survey was conducted as well (Figure 6). 

This data was used to determine the pseudo-vorticity 


Circulation of the region using the C-Vector concept. 


A. DATA COLLECTION 
ie) CLD 
All CTD stations with their positions, weather 
conditions and maximum depth (m) are shown in Table II. 
Stations 1, 20, 21, and 40 were not included because they were 


located outside the grid (Figure 9) which was used for the 


Z3 


Table II.CTD STATION LOCATION, WEATHER AND MAX.DEPTH 





Stat. Lat. Sas Se “Wind | Air Max. 
No. Deg. Deg Min Dir. (m/s) Temp. Depth 
N (C) (m) 
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calculation of the C-Vector. The CTD casts were made with a 
Neil Brown Mark III-B CTD attached to a General Oceanic 
rosette sampler for in situ water sampling. This model CTD is 
considered to be accurate to within + 0.005 mmho, + 0.0059 ¢, 
emiad t+ 6.5 db. Conductivity, temperature and pressure 
resolutions are 0.001 mmho, 0.0005° C and 0.1 db respectively. 
The CTD sensors were calibrated prior to the cruise. 
Coefficients calculated for these calibrations were applied to 
the raw data. Salinities were further calibrated by adjusting 
the calculated salinity to the salinity of water samples 
collected during the cruise. Further details of instrument 
calibrations can be found in the NPS data report for this 
Cruise (Jessen et al., 1992). 
2. ADCP 
Continuous relative current velocity data was gathered 
aboard the R/V Point Sur using RD Instruments Vessel Mounted 
Acoustic Doppler Current Profile (VM-ADCP), equipped with a 
four beam JANUS array, operating at a frequency of 150 khz. 
3. Wind data 
During the cruise, an R. M. Young anemometer on the 
R/V Point Sur obtained 30 second averages of wind speed and 
direction (Table II). 
4. Satellite 
The AVHRR Instrument on NOAA-11 satellite measured 


emitted radiation of the study area in five wavelength bands 


aS 


(three thermal Infra Red (IR) Bands, one Near Infra Red (NIR) 


and one Visible Band). 


B. DATA PROCESSING 
Ta cLD 

The raw CTD data was processed on a PC using an EG&G 
Marine Instrument software package called CTDPOST. This 
program flags bad data values using a first difference 
criterion and allows the user to examine and interpolate 
across bad data if necessary. After the correction of the 
data, a dual time lag filter was applied to remove time lag 
spikes, and the data was averaged to 2 db. 

2. ADCP 

The conversion of recorded ADCP relative velocity data 
into absolute velocity data took place on board ship. This 
process depends significantly on the accuracy of the ship’s 
navigation instruments. The ship’s position data was obtained 
from GPS and was also recorded on computer disks for later 
processing. 

Variability of the ship’s speed, direction, and data 
collection interval, play significant roles in the ADCP 
accuracy and reliability. More details concerning ADCP data 
processing techniques can be found in Kosro (1985). According 


to Kosro, accuracies of 4 - 5 cm/sec in the athwartship 
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component and 2 - 4 cm/sec in the fore-aft component are to be 
expected. 
3. Wind data 

The wind data, collected underway, was stored on the 
hard disk of an HP310 computer. This data was transferred to 
5.25 inch diskettes upon return to port and computer processed 
to obtain hourly averaged wind vectors (Figure 4). 

4. Satellite data 

The Satellite data was first calibrated, transferred 
to earth coordinates and temperature calculated pixel by 
pixel. All images were coregistered to a Similar coordinate 
map of the study area. This work took place at the Naval 
Postgraduate School IDEA Lab. 

5. Geostrophic velocities 

All densities and geostrophic velocities were obtained 
from the calibrated CTD data uSing computer programs developed 
in the NPS oceanography department based on the equation of 
State of seawater (Fofonoff, 1985). 

The u and v components are the across shore and along 
shore components of the geostrophic velocities. These 
components at a particular station were calculated using 
Stations on either Side of it. For example, utilizing 
Stations 26 and 24, the v-component of geostrophic velocity of 
Station 25 can be determined, and from stations 16 and 36, the 


u-component of station 25 (Figure 9). 
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Figure 9. Distribution of ClD@stausenuce 


This does not work for the external stations. In this 
case the u and v components of the external stations were 
calculated using the external stations and the stations close 
EO Pehem. For example, for the determination of the u- 
component of geostrophic velocity at station 3, stations 3 and 


18 were used. 
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6. Reference level or Level of No Motion (LNM) 

The calculation of geostrophic velocities depends on 
the right choice of a reference level or Level of No Motion 
(LNM) . 

The hydrographic and ADCP data from the NPS Farallones 
Shelf and Slope study on May 16 to May 21, 1992 were examined 
carefully. It was determined from the plots of the vertical 
sections of 5 km averaged, across and along-transect ADCP 
velocities (for sections A, B, C, D and E) that the existence 
of a common horizontal reference level or LNM was uncertain. 
It was obvious that the reference level could not be 
horizontal over a great distance. 

In this study a large number of stations were located 
in shallow waters so the area was divided into two parts 
(Figure 10). A shallow part which includes stations with 
depths less than 1,000 meters (over the continental shelf), 
and the deep part which includes stations with depths more 
than 1,000 meters (over the continental slope). It was 
assumed that the reference surface was horizontal within the 
deep waters. 

According to Soule (1939) the water density field has 
to be analyzed first and the horizontal reference surface has 
to be placed in the layer where density is least disturbed. 

In analyzing the distribution of density differences 


between neighboring stations the depths were found where 
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Figure 10. A 3-D view of the stations over the shallow an 
deep waters. 


density differences changes sign or became zero. These depths 


represent the levels of minimum current shear. 


Figure 11, which illustrates the plots of the densi, 


differences between neighboring stations in the shallow part 


of Section D, shows that there is no common depth where the 
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Figure 11. Density differences between pairs of stations of 
Section D. 


density differences are zero. The choice of the LNM in 
shallow waters is therefore uncertain, and the reference layer 
was estimated to be the lower level of measurements between 
neighboring stations. 

The analysis of density difference distribution for 
deep stations offers better results. Figure 12 includes all 
the pairs of deep stations. A common depth of "zero" density 
differences is found at 1,400 meters. This level (1,400 m) 


was used as a reference level for deep waters. 
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There is another method of determining the LNM similar 
to the above process (Mamaev, 1955). The difference is that 
the curves for the vertical distribution of differences of 
specific volume at two neighboring stations are analyzed, 
instead of the differences of density. 

7. Computation of C-vector 
decal ciuijerommons C-vector (C,, €,, Cz) from the 
geostrophic currents. 
The computation was done using equations (16) - 
(18), where: 
® U, and V, can be calculated for all stations. 
@® Coriolis parameter was assumed constant, f=2Qcos¢ (where 
Oe eo oe eee sean an average latitude ¢ = 37° 30’ 
North 1s used) 


@ YY and Y~, the turbulent momentum fluxes due to surface 
wind stress (17*,7TY), are given by equations (6) 


Taking the average density of all stations it was 
@ecermined that an average ocean density of the region was p= 
1.026 gr/cm?. 


WonGact sess) as Obtained from 


t%=p ,C,| Vu, ’ v =p, ©-| Vale va) 
m@ere p, = 1.293 103 gr/cm? (iewadenotey)yec, = 1.3 10°73 
(drag coefficient), |V,| is the hourly averaged wind speed at 


Sr 


height 10 m and u,, v, are the x and y - components Givi 
speed respectively. 

The mixed layer depth (MLD) of the study area is 
needed to determine to what depth the wind mixing influences 
the ocean stratification. Observing all MLD from temperature 
plots during this study, it was determined that an average MLD 
(of all stations) is 15 m. This means that wind stress mixes 
the waters down to 15 m. 

A linear reduction of wind stress with respect to z 


was assumed (Figure 13) for the purposes of C-Vector 


Calculation. 





Figure 13. Linear reduction of wind stress with respect to 
depth. 


b. Computation of C-vector |(Gage C,) due to Weae 


y! 
total flow (geostrophic and ageostrophic). 
This computation depends on the equations (22) - 
(24) where C,, C, and C, and geostrophic velociticomie. 


calculated as described above. 
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Having C,, and C, it is then easy to compute W from 
equariGmm 29) Contours of C-Vector for ageostrophic flow, C- 
Vector for the total flow and W are analyzed and compared in 
the next chapter. These three parameters are the main tools 
for the identification of the ageostrophic and the total 


(geostrophic and ageostrophic) pseudo-vorticity fields. 
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IV ANALYSIS OF PSEUDO-VORTICITY FIELDS 


A. DESCRIPTION OF THE ANALYSIS 
This chapter presents a detailed description of the 
pseudo-vorticity fields for the ageostrophic cCircWiaeaa7 
(C,,C,,C3) in the Gulf of the Farallones. The results Geamae 
description enable us to proceed further in the analysis of 
the pseudo-vorticity fields for the total geostrophic and 
ageostrophic circulation (C,,Cy,C,). Then the vorticity of C- 
Vector (WV) will be used to indicate the upward and downward 
motions in the area, and it will be compared to a satellite 
picture of the same area and the vertical ADCP velocities. 
1. Pseudo-vorticity fields (C1 ,C5,C3) for the 
ageostrophic circulation 
a. Analysis of the x-component of ageostrophic 
pseudo-vorticity field 
C, is the ageostrophic pseudo-vorticity component 
of C-Vector in the x-direction or cross-shelf direction. This 
direction approximately indicates the cross polar vertical 
circulation pattern. The values for Ga ee have been plotted 
on the y-direction or along-shore direction. 
From Figure 9, notice that there are nine columns 
on the y-direction. Columns 1, 2, 3, 4 and 5 corresponaiias 


shallow water and are going to be analyzed first. 
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Puta ole > ana l6esnow the, plots of the first 
five columns. The arrows on the plots indicate the direction 
of rotation of the ageostrophic pseudo-vorticities and not the 
Girections of currents. The cells with positive (negative) 
values of the ageostrophic pseudo-vorticity are turning 
counter-clockwise (clockwise) or cyclonic (anti-cyclonic). The 
Girection along the positive (negative) axis is upward 
(downward) out of the page (into the page). 

From these five plots it can seen that the 
ageostrophic pseudo-vorticities are increasing going offshore. 
At column 1, the maximum values of ewer are about +/- 0.1 and 
at column 5 the values have increased to 1.3. The first two 
columns have the negative ageostrophic pseudo-vorticity cells 
north and the positive cells south. The opposite occurs at 
columns 3 and 4 and at column 5 the positive ageostrophic 
pseudo-vorticity predominates from South to North. 

At column 5 the strong south positive ageostrophic 
pseudo-vorticity cell coincides with the area of large 
temperature and salinity gradients of Section A (Figure 6). 
Also, the strong middle positive ageostrophic pseudo-vorticity 
cell coincides with the significant salinity gradient of 
Section D (Figure 6). Column 5 has the largest positive value 
of the ageostrophic pseudo-vorticity fields. 

Pqunes Loyal 7scanad .28 show the four plots .of 


columns 6, 7, 8 and 9 which lay over the deep waters. Column 


ay) 


7 shows a strong positive ageostrophic pseudo-vorticity field 
that dominates the area. Columns 6 and 8 have their negative 
ageostrophic pseudo-vorticity cells on the north side and 
Column 6 has one positive cell south. Column 8, however, has 
a positive cell followed by another negative cell south of the 
North weesitive celle Column 9 has a significant negative 
north ageostrophic pseudo-vorticity field. Its position 
coincides with the large temperature and salinity gradient of 
the west side of Section E (Figure 6). 

The significant temperature and salinity gradients 
are probably due to subarctic water underlying a surface lens 
of fresh and warm water (Collins, 1992). Figure 19 presents 
a better picture of the ageostrophic pseudo-vorticity fields 
through a 3-D plot of columns 9,7 and 5. The red lines are the 
solid positive values and the blue lines are the dashed 
negative values. 

The majority of the ageostrophic pseudo-vorticity 
fields are located at the upper 200 m. Waters below 200 m 
were found to have very weak ageostrophic circulation, so the 
significant ageostrophic circulation is close to the surface. 

b. Analysis of the y-component of the ageostrophic 

Ppseudo-vorticity fields 

The values for C,/f? have been plotted at the five 
Sections A, B, C, D, and E (Figures 20, 21 and 22). Section 


A shows a relatively weak ageostrophic pseudo-vorticity field. 
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GiewemeOrSsso=® Secer1onseB and GC include four cells. Moving off 
shore alternating values of negative and positive cells are 
observed. Sections D and E have the highest values of the 
non-dimensional component C,/f*. The west (left) side of 
Section E has the maximum negative value of the ageostrophic 
pseudo-vorticity fields. This location coincides with the 
large off - shore temperature and salinity gradient of section 
E (Figure 6). The values of Crip below 200 m are very small 
as can be seen from the plot of Section C (Figure 21). 

The magnitudes of Ga /ere (maximum value -7) are 
larger than the magnitudes Ciyase (maximum value 1.3). This 
probably happens because the average distance between the five 
sections (dy = 20 km) is bigger than the average distance 
between the nine columns (dx = 8 km), and the derivatives with 
respect to y are smaller than the derivatives with respect to 
oe. 

Figure 23 offers a three dimensional view of the 
ageostrophic pseudo-vorticity cells of Sections A, C and E. 
All the sections show that the ageostrophic circulation is 
located close to the surface and close to the locations of the 
large temperature and salinity gradients. 

c. Analysis of the z-component of ageostrophic 
pseudo-vorticity field 

As expected the values of vertical component yt 


are smaller than the values of the horizontal components 
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because the vertical ageostrophic vorticity is much smaller 
than the vertical geostrophte Vvorereme. 

The upper plot of Figure 24 shows Cae at the 
surface, without the wind affect. Notice that there are tae 
negative ageostrophic pseudo-vorticity cells and one positive. 
The large values of Caer are concentrated at the north-west 
Side of the area where the large temperature and salinity 
gradients were observed before. The lower plot of Figure 24 
shows the contours of Gyre” including the wind. The majority 
of the values of Cie: are negative at surface. This probably 
happens because the area with the negative values coincide 
with the area of strong winds. A comparison of the plot of 
Cry ban at the surface with wind and without wind shows that 
without wind there is one cell with relative large positive 
values of Ca) fe on the North-West side. 

The plot which includes the wind shows that the 
negative values predominate in the region, and the positive 
cell that existed without the wind has become smaller and 
weaker (the magnitude is smaller). For Cy) at 50 m (Figure 
25), the three negative cells and the one positive cell exist 
at the same positions but with smaller values. 

At 200 m (Figure 25), the three negative cells 
still exist almost at the same location but with smaller 
magnitudes. The previously formed positive cell is very weak, 


and a new positive cell appears at the south-east corner of 
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the area. Figure 26 shows a three-dimensional picture of the 
three plots of Cyt? at the surface and the depths of 50 m and 
2100! tm. 

Figure 27 shows a 3-D picture of C-Vectors of 
columns 9, 7 and 5. The large negative values of C-Vector 
(for ageostrophic circulation) were found at the west side of 
Section E where the large temperature and salinity gradients 
were observed. 

2. Pseudo-vorticity fields for the total circulation 
(C,,Cy,C,) ageostrophic and geostrophic 
a. Analysis of the x-component of C-Vector (C,) of the 
pseudo-vorticity field 

The columns 1, 2, 3, 4 and 5 are located over the 
shallow waters. At these columns the negative pseudo- 
vorticity fields are dominant (Figures 28, 29 and 30). The 
strongest cells are located south and north of column 5, and 
the south one coincides with the large temperature and 
Salinity gradients at the middle of Section A. 

At the deep waters, the negative cells of pseudo- 
vorticity are dominant also (Figures 30, 31 and 32). At the 
north part of column 8 and 9, the maximum values of C,/f are 
observed (C,/f = -35). In the same area, maximum values of 
ae f° and Cae and Sy eae were seen before. This is the area 
with maximum temperature and salinity gradients. The arrows 


show an upward flow at column 9. At column 8, there is an 
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upward motion south and north and at the middle downward flow. 
The arrows indicate the upward and downward motion of all 
columns. 

Figure 34 is a three-dimensional plot of Cy aaa 
columns 9, 7 and 5 where the upward-downwards motion can be 
seen more clearly in three-dimensional space. 

b. Analysis of the y-component of C-Vector (Cy) of the 
pseudo-vorticity field 

Positive pseudo-vorticity cells of Cy dominate the 
five sections (Figures 34, 35 and 36). Notice that a negative 
pseudo-vorticity cell exists at the west side of all Sections. 
The arrows show an upward motion off shore of all Sections 
followed by a downward motion (going east). At Sections A, B, 
and C, there is an indication of upward motion at the east 
Side of these sections. Sections D and E have the strongest 
pseudo-vorticity cells at their west side (maximum Cy /f = 13 
ate Section Db), This area contains the large horizontal 
temperature and salinity gradients. Figure 37 contains a 
better picture (3-D) of the pseudo-vorticity fields of 
Sections A,C and E. 

c. Plots of ¥Y function 

Figure 38 shows W at the surface. The large 

positive values of VW mark an upward motion at the west side 


of the region. In the same figure, there is an outer area 
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west of the Farallones Islands where there are positive values 
of W and these are probably upwelling waters. 

Between these two areas, there iS an area with 
negative values of VW (downward motion). This area lies over 
the continental slope. A weak upward motion (positive values 
of W) appears at the east side of the region and south-east 
of Farallones Islands. At 20 m, the picture is similar to the 
surface. 

At 50 m (Figure 39), the locations of upward and 
downward regions are almost the same as those seen at the 
surface and 20 m. At 200 m, the situation starts to invert, 
and the regions of upward motions become regions of downward 
motions. This means that the upward-downward motions do not 
extend as deep as 200 meters. 

Figure 40 is a 3-D plot of W at surface, 50 meters 
and 200 meters. iimeeits wecCloOmmoLcture gthe solid positive 
values are red (upward motion) and the dash negative values 


are blue (downward motion). 


B. VERIFICATION OF C-VECTOR VORTEX LINES 
1. Direct verification of C-Vector vortex lines by WV 
function 
An important verification of the upward and downward 
water motions at Gulf of Farallones Islands induced by C- 


Vector can be achieved by the W function. The positive values 
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of W give the upward motion (upwelling) and the negative 
values the downward motion (downwelling). 

From examining the plotsiior |) ese Cy it was 
determined that when C, and C, vortex lines indicate an upward 
(downward) motion, they were consistent with the positive 
(negative) W lines. 

2. Indirect verification of C-Vector vortex lines by SST 
images 

Figure 41 1S an AVHRR image of the area for the day 
prior to the cruise (15 May 1991). It suggests that a plume 
of cold water (8° C) extends offshore of Point Reyes and over 
the continental slope. This area coincides with positive 
values of W (upward motions) south-west of Farallones Islands 
(Figure 38). 

South-east of the continental slope and north-west of 
the Pioneer Seamount, there is an area of relatively warm 
water (10° - 11° cc). At this area, W has negative values 
(downward motion). At the south-west side of Pioneer 
Seamount, there are relatively cold waters (9° - 10°) andi 
the same area, W has positive values. Therefore, the SST 
images provide useful information concerning possible vertical 


velocity at the ocean surface and are a helpful indirect 


VeErlEleGaelon of v. 


44 


3. Comparison of vertical ADCP velocities with Vv 
function 

Figures 42 and 43 shows the vertical ADCP velocities 
at 20 m and 50 m at the upper part of the figures, and WV at 
20 mand 50 m at the lower part of these Figures. Notice that 
at 20 m depth and at the south-east of the region, both WV and 
vertical ADCP velocities are consistent because they show 
negative values. The W has negative values at the center of 
the region, but the vertical ADCP velocities are negative at 
the south-west side of the area. These two regions of 
negative values coincide partially. The rest of the areas 
have positive values. 

At the plots of 50 m, the vertical ADCP velocities 
indicate negative values at the same areas as at 20 m. The ¥ 
shows negative values at the east side of the region and at 
the middle and south region. The areas of negative values of 
W and negative vertical ADCP velocities partially coincide and 
partially differ. 

There is no perfect match. This happens probably 
because the vertical ADCP data are also affected by tides, 
waves or internal motions. Also, the calculation of VY 
contains some small errors like the round-off error, 


truncation error and machine error. 
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Figure 24. C,/f* values at surface without wind (upper pice) 


and with wind (lower plot). 
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Figure 43. Vertical ADCP velocities at 50 meters 


and ¥ function at 50 meters 


(Mower plot) . 


(upper plot) 


V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

The purpose of this study was to demonstrate the C-Vector 
method and the potential value of this method in the 
calculation of the ageostrophic circulation using CTD and wind 
data. Sometimes the ageostrophic flow is significant and the 
purpose of this study 1s to emphasize the importance of 
ageostrophic circulation in waters off California as inferred 
by C-Vector using CTD and wind data. 

Spring is an upwelling season for waters off northern 
California. The SST satellite image of the area shows two 
filaments of cold water. One extends west of the Farallones 
islands and the other stretches offshore of the 1,000 m 
1sobath. The cold filaments are associated with positive 
values of the vorticity of C-Vector W (equation (29)). 

The main upwelling extends to a depth of 200 m, and the 
strong C-Vector cells extend from the surface to about 200 m. 
Cells below 200 m are observed to be very weak. This is in 
phase with the Sverdrup assumption that the upwelling water 
rises from moderate depths only, probably less than 200 m 
(Sverdrup and Fleming, 1941), and the phenomenon represents 


overturning of the upper layer. 
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Some large values of C,, Cz, and C,, Cy, (horizontal 
components of C-Vector for the ageostrophic and the total 
Circulation) are found over the continental slope. This 
probably occurs because of the presence of the continental 
slope and is associated with oceanic processes which occur 
there. 

A shortcoming in the calculations of C-Vector could be the 
uncertainty associated with the determination of a reference 
level. Density differences between stations demonstrate that 
1,400 m is the main representative LNM for the deep waters. 
While the depths where the difference in density between 
Stations 1S zero are not always 1,400 m, there is a good 
indication that at this level currents are very small. 

The Beta-spiral method (Schott and Stommel, 1978) 


evercomes the erroneous "level of no motion" and allows the 


"level" to vary realistically in (x,y). 


B. CONCLUSIONS 

The C-Vector is a good method for calculating the three- 
dimensional pseudo-vorticity fields in the ocean using CTD 
measurements and wind stress data. 

This technique (with the help of W) provides some useful 
information of the upward and downward motions in the ocean. 

If the pseudo-vorticity equations can be solved, the 


three-dimensional circulation in the ocean can be obtained. 
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C. RECOMMENDATIONS 

Data from May 1991 gave a picture of the ut1lityeeneee 
Vector method and the pseudo-vorticity fields at the Gulf of 
Farallones Islands. 

These data are complicated by such local phenomena as 
upwelling, tides, inertial motions and winds. T teas 
recommended that the results of this study should be compared 
Hebe) ee 

@® The pseudo-vorticity fields using data from other sources, 
for example from a current meter array covering the area 


for verification 


® Results of other studies and calculations of C-Vector teem 
the same area. 


In the future it is also recommended to use the 6-spiral 
method in the C-Vector computation. The advantage of the 6- 
Spiral method is that it overcomes the problems of the "level 
of no motion" and for this reason it is expected to improve 
the C-Vector method. 

Another way to avoid the uncertainty of the level of no 
motion or the reference level could be the use of satellite 
altimetry data. This data can be used for the determination 
of the geostrophic currents at the ocean surface. Then, an 
integration of the thermal wind equation from surface to a 
certain depth z can provide the geostrophic velocities at 


different depths. 
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